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SUMMARY 



Simple curves are presented by which the basic 
pressure-drop characteristics of unhealed tubular 
radiators can be corrected to operating conditions in 
which the radiator is heated and in which the Kach number 
of the tube flow is of appreciable magnitude. The only 
data required for the use of the curves are the radiator 
dimensions, the rate of heat input, the pressure and 
temperature ahead of the radiator, and the rate of mass 
flow of air through the radiator. 

The accuracy of the curves for predicting the 
compressibility effects for unheated radiators is 
confirmed by comparison with test data obtained from two 
.independent sources. An example of the use of the curves 
for a typical oil-cooler installation is given. 



INTRODUCTION 



The pressure-drop coefficient for a cold tubular 
radiator, as usually determined from tests with low 
airspeed in the tubes, is known to be subject to 
correction if the radiator is operated under conditions 
in which an appreciable density decrease occurs as the 
air passes along the tube. This decrease in density may 
be caused by the increase in absolute temperature due to 
the addition of heat, by the reduction in absolute static 
pressure due to the pressure drop within the tube, or by 
both. The density decrease due to the addition of heat 
is important under all operating conditions, but the 
density decrease due to the static-pressure reduction 
becomes important only when the pressure drop is appreci- 
able in comparison with the absolute static pressure - 
that is, when the Mach number of the tube flow becomes 
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of appreciable magnitude (reference 1). In high-speed 
airplanes, neglect of the Mach number effect will result 
In sizable underestimation of the pressure drop required 
to Induce the necessary cooling-air flow, 

A number of investigators have demonstrated that 
comoressibili ty has no appreciable effect on the heat- 
transfer coefficient. The temperature difference on 
which the heat-transfer calculations at high Mach numbers 
must be based, however, is the difference between the 
wall temperature and the stagnation temperature of the 
cooling air rather than the true temperature of the 
cooling air. (See reference 2.) Test data on the heat- 
transfer characteristics obtained at any test Mach number 
may thus be used to determine the mass flow required for 
cooling for the design condition. In the present report 
the mass flow is assumed to have been so determined 
and the report is concerned only with the determination 
of the pressure drop necessary to Induce the required 
mass flow. 

The purpose of the present paper is to present 
curves from which the basic pressure-drop characteristics 
of cold radiators can be corrected to operating conditions 
in which the radiator is heated and the Mach number of 
the tube flow is of appreciable magnitude. The only 
requirement for the use of the curves is knowledge of 
the radiator dimensions, the rate at which heat is to be 
dissipated, the pressure and temperature ahead of the 
radiator, and the rate of mass flow of air through the 
radiator. The curves are constructed for a Mach number 
range from zero to the maximum attainable Mach number 
corresponding to existence of sonic velocity at the 
exits of the tubes. 

A knowledge of the theory by which the curves are 
derived is not essential to their use. ?or the sake of 
completeness, however, this theory is briefly reviewed 
herein. For the convenience of the reader, the report 
is presented under three main headings. The theory is 
found in the section entitled "Analysis. " The validity 
of the theory and the accuracy of the curves is confirmed 
by comparison with experimental data obtained from 
references 3 and 1+ in the section entitled "Comparison 
of Theoretical Results with Experimental Data." The use 
of the curves is illustrated by a typical example in the 
section entitled "Example of Use of Curves." 
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SYMBOLS 



A cros3-rectional area of duct or radiator tube, 

square feet 

a velocity of sound in air, feet per second 

Cf skin-friction coefficient (Vf/^rz^) 

Cr) f mean skin-friction drag coefficient of radiator 

tube (Df/q r2 Ar 2 ^ 

c- 0 specific heat at constant pressure (for air, 

0.21a Btu/lb/°F) 

Df drag force due to skin friction in a radiator 

tube 

d radiator-tube diameter, feet 

g acceleration of gravity, feet per second per 

second 

H beat added in radiator, Btu per second 

h total pressure, pounds per square foot 

Ah total-pressure loss, pounds per square foot 

L radiator-tube length, feet 

M Mach number (v/a) 

m mass-flow rate, slugs per second 

m cr mass-flow rate at whic^h sonic velocity is attained 

at tube entrance I - = 0 ) 

\d / 

p static pressure, pounds per square foot 

Ap static-pressure decrease, pounds per square foot 

q dynamic pressure, pounds per square foot (ipv/ 
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R Reynolds number or 

S area of inside surface of radiator tube, square 

feet 

T air temperature, °F absolute 

AT air temperature rise, °P 

v velocity in radiator tube or in duct, feet per 

second 

p density, slugs per cubic foot 

F c compressibility factor (l + |^ + |j~ + • . 

\i viscosity of air, pound-seconds per square foot 

Subscripts (see fig. 1): 

2 station in duct ahead of radiator 

J station in duct behind radiator 

f friction component 

i low speed, incompressible-flow condition 

V2 within radiator at tube entrances 

rz within radiator at tube exits 



ANALYSTS 



Derivation of Equation Relating Pressure Drops for 

Compressible and Incompressible Plow 

For the purpose of the present analysis the pressure 
drop across the radiator is taken as the difference 
between the total pressure at a station just ahead of 
the front face of the radiator and the static pressure 
at the tube exits. A diagram of the flow system across 
the radiator is shown in figure 1. The station desig- 
nations used in this figure were taken from reference 1. 
The pressure dron across the radiator for this system is 
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defined by the equation 



Ap - h^ - p r 



It will be shown subsequently herein that the results of 
the analysis ozn bo aoplied to other definitions of the 
pre 5 a are d^o-, >var,h as p 2 - Pj> ^2 " p r 2| > or h 2 " h 3 

The pressure drop can be conveniently divided into 
two cai.r."o:\^n-£; ; the pressure drop due to acceleration 
of the iur into the tube entrance and the pressure drop 
within the tube duo to friction and momentum change. 
From equation 13 of reference 1, Ap may be expressed as 
the sum of these two components: 



This equation can be evaluated with the aid of data given 
in references 1 and 5, but the orocess is not simple. The 
problem can be simplified by assuming that the pressure 
drop for cold Incompressible-flow conditions, herein 
designated is known from basic test data for the 

radiator. The value of Ap for operating conditions 
can then be obtained if the ratio Ap/Ap^ is evaluated. 

Vrom equation (1), for the incompressible unheated 

/ P ro \ 

condition ( F c - 1, — - = 1.0 j, the calculated value of 



Ap .- ]Qr2 




(i) 




(2) 



More specifically, since 




in incom- 



pressible flow and if 




is replaced by 



CD fl , 




CONFIDENTIAL 



6 



CONFIDENTIAL NACA ACR No, lij.123 



An approximate application of equation (2a) that is 
in general use for estimating radiator pressure drops for 
conditions other than the test condition is 

P2 test 



Ap ~ Ap test— 



The "test 11 pressure drop is measured with the mass flow m 
required for the design heat dissipation. This formula 
obviously neglects both Reynolds number effects (changes 
in Cj^. ) and compressibility effects, both of which are 

too large to be neglected in many present-day applications, 

It is convenient at this point to list the equations 
for evaluating Ap* for all usable definitions of the 
radiator pressure drop. Fro.n equation (2) or (2a), 



h 2 - . r 



3l - 1 + Cj>. (2b) 




Similarly, since p 2 = h 2 - q 2 ar i d <l2 = \p£ 



ho ~ Pz\ 

In rd r be e aluatt " v y an assumpti n must oe 



made as to the location of station J. It is assumed that 
this station is far enough behind the radiator for the 
velocity distribution across the duct to be uniform. Then, 
from equation 12 of reference 1, for incompressible flow, 



( 



v *r 2 /. " >3 " \*3 / | 
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whence , 



1 + C D „ - ( — -2!^ 



/A \ ^ 

/ A r j 



As 



(2d) 



The solution for incompressible flow. 

h 2 - h 5 = d 2 - "5 + q 2 " q $ 

J q ^2 \ q - V-2 J ' V l 3 / 

/h 2 - h,\ 

is used to evaluate — by use of equation (2d) as 

\ lr 2 / 



h 2 " h 3 \ , „ * r 2 / Ar 2 



2 



\ 



, = 1 + C D ~ - 2t-^ + ( — --) (2e) 



As discussed in reference 1, the quantity Df . Mr^n 

or Cj)^ Is a drag coefficient. The value of this basic 
1 i 

friction factor depends on the skin-friction coefficient 
and on the tube length-diameter ratio. From the definition 

Of Cf , 

The skin-friction coefficient c f for a given tube 

deoends only on the Reynolds number (references 3 and 1+) . 
The addition of heat within a tube has little effect on 
the local Reynolds number because the quantity pv remains 
constant along the tube and, in usual cases, there is 
little change in viscosity. Figure 2, which is based on 
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the relations presented In reference 5> allows Cf 

and Ctv-. to be obtained for any tube L/d ratio for a 

wide range of Reynolds number or mass-flow rate. As 
indicated in figure 2, the Reynolds number may be 
evaluated directly from the mass-flow rate, the area, and 
the viscosity. For usual applications the viscosity at 
station 2 may be used instead of the viscosity at 
station ?z with negligible error in the resulting 
pressure -drop evaluation . 

The ratio of equation (1) to equation (2) can be 
expressed as 




Prom reference 1, if q£ ia sma11 iri comparison with q r £, 



p r 2 



= ( 1 + 0.2lv' r9 2 Y° 



Therefore, 



An = _( 
Ap< 



1 + 0 • 2Mpp^ N \ 



i + c Dfi 



crp 



D f 



CD fi + ^ 



'Pr 2 



i) 



(3) 



Evaluation of Equation (3) 

The method of evaluation of the terms in equation (3) 
for any combination of the known conditions ahead of the 
radiator, the mass-flow rate, and the heat-inout rate will 
now be discussed. 
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The I,!ach number at the tube entrance U r ^ 9 which 

aopears in equation (3), depends only on the known mass- 
flow rate, the known temoerature and pressure at 
station 2, and the leak area through the radiator. From 
figure 3, can be obtained directly for any 

combination of these basic conditions. 

The terms Co f > Cr )f , and P^/Pr^ in equation (3) 

are interdependent. The friction drag coefficient C Df 

for the actual operating condition is greater than CDf., 

except when the tube flow is cooled appreciably, as in 
the case of charge-air flow through intercooler s . The 
ratio CDf/CDf, is dependent on the values of M r ^ 

and Pr 2 /Prj« These values govern the increase in dynamic 

pressure along the tube and hence determine the drag 
force due to skin friction Df under actual operating 
conditions. In reference 1, an approximate empirical 
relation was assumed, as follows: 




This relation has been found to be inaccurate except 
for values of P^/Pr^ n ^ar unity. An exact evaluation 

of CD f /CD f has been made for unheated tubes with 

the aid of reference 6. Equation 12 of reference 6 
gives the following relation between CDf., Kr?* 

and Pr 2 /Pr 5 , obtained from a solution of the differential 
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equation of the unheated flow through a constant-area 
tube • 




Prom equation (10) of reference 1 the corresponding 
relation for Co f > based on considerations of continuity 

and of momentum and energy changes between the entrance 
and exit of the tube, is 




The ratio of equation (5) to equation (J4.) has been 
evaluated for the complete range of density-ratio values 
corresponding to a number of constant values of 

and the results are plotted in figure I4.. 

The procedure then used in the evaluation of Cr^ 
and P^/^rj with the aid of figure I4., for given values 
of C Df> and M r ^, was as follows? 

(1) On the assumption that Cj^ ~ CDf^* equation 5 
was used to obtain an approximate value of p r ^/p r ^. 
Figure 5(a) of reference 1 is a solution of equation 5. 

CONFIDENTIAL 



NACA ACR No, ll|l23 



CONFIDENTIAL 



11 



(2) With this value of Prz/Pry c D f / c D f 



was 



determined from figure l± and a second approximation 
for Cj) f was thus obtained. 

(3) This more exact value of Cj^ was used and 
the procedure was repeated to obtain a third approximation 
for C Df and Pr 2 / p rj • 

The values determined from the third approximation 
were generally found to have the desired accuracy. 



For unheated tubes all the quantities involved in 
the pressure-drop ratio (equation (3)) have been shown 
to be functions of only the two fundamental variables C 



and ^- r2 - The variable depends on the tube 

dimensions and Reynolds number, and the variable 

depends on the mass flow and state of the air ahead of 
the radiator. Both variables can be simply evaluated 
for a given installation, as previously explained, for 
the particular design conditions under consideration. 
Equation J has been evaluated for a range of these 
oarameters and the results are olotted in figure 5- 
Figure 5(a) shows the pressure-drop ratio plotted against 

for unheated radiator tubes of various CDf^ values. 

When heat is added to the tubes, the density ratio 
p r ^yPr^ is increased and the drag coefficient Cj>f is 

correspondingly increased for a given mass flow. In the 
evaluation of CD f for the heated condition, the 

ratio Q'DffG'Df from figure ii. was .assumed to be the 
same, at a given p r ^/p r ^ and M r ^, as for the unheated 

condition for which figure i| is exact. This assumption 
is believed justified for practical purposes. If heat 
were added uniformly at low speed, linear reduction in 
density due to heating would occur along the tube. For 
a given value of P^/Pr^ > the ratio c Df/ c Df£ would be 

greater for a heated tube than for an unheated tube 
because of the nonlinearity of the density variation for 
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the unheated tube at the higher Mach number at which the 
unheated tube would have to operate to give the same 
p r2 /p r ^ value as the heated tube. Results based on the 

assumption that the drag-coefficient ratio is the same 

for heated or unheated flow will thus tend to underestimate 

the required pressure drop when heat is added. The 

probable error will be seen to be small for usual 

operating conditions because the oressure-drop ratio Ap/Ar^ 

is not critically dependent on the value of CDf for 
these conditions. The tube Mach number for these conditions 
is generally well below the critical, and the density 
variation along unheated tubes then approaches linearity. 
The assumption made in calculating CDf for the heated 
condition would therefore not be expected to result in 
appreciable error when the entrance Mach number is well 
below the critical. 

The relationship among CDf, Mr 2 > and Pr2 /Prj 
for heated tubes is given in figure 5 of reference 1 for 
constant values of the heat-Input parameter c g^T r2 ' 

This figure was used with figure It as described for the 
unheated condition to obtain corresponding values of CDf ± > 

C Df , and Pr 2 /Pr>3 ?or constant values of Cpg S Tr2 - 

Values of Ap/Api were then computed from equation 5 and 
are shown in figures 5(b) to 5(g)- Sach part of figure 5 
corresponds to a constant value of c g^T^ ' As dlscussed 
in reference 1, for the incompressible condition, 

H _ AT 
c P gmTr2 

Interpolation will be required to determine Ap/iPi for 

values of the heat-input parameter intermediate to those 
for which the curves were derived. The heat-input 
parameter may be evaluated from the following relation 
derived on the assumption that v 2 is small in comparison 
with vr^: 

H _ /' H \ (, n 2\ (6) 
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The results shown in figure 5 indicate that appreciable 
compressibility effects on the radiator pressure drop 
exist even at relatively low entrance Kach numbers, 
particularly for tubes of large values of L./d (large 
values of C Df ^ ) and for installations in which 

high rates of heating are used. The limiting values 
shown on the curves correspond to the attainment of sonic 
velocity at the tube exits. No increase in mass flow 
through the tubes can be effected by lowering the exit 
pressure below the Uniting value corresponding to sonic 
exit velocity. It is interesting to note that this 
limiting condition occurs at M r ^ = 0.1+8 for an unheated 

radiator tube of Ccf^ = 1.2 (fig. 5(a)), which corresponds 

to a 0,196-inch diameter tube with a typical value of L/d 
of about 60 for standard conditions ahead of the radi- 
ator ( fig. 2) . 

The curves of figures 5(h) to 5(g) include the 
effects of both heating and Mach number; that is, the 
figures as given show the net effect of compressibility 
(density change) on the pressure drop. If separation of 
the two effects is desired, the Mach number effect alone 
can be determined from figure 5( a ) f° r comparison with 
the net effect for the heated condition. 



General Application of Equation (3) 

Test data Tor radiators are frequently given in 
terms of - P j y V? - p r ^, or hp - hj, rather than 

in terms of Y12 - Pr^> which was used in the present 

analysis. The absolute values of the pressure drop vary, 
depending on the definition, by 5 to 15 percent in usual 
cases. The ratio Ap/Zio^, however, can readily be shown 
to be essentially the same for all the definitions of the 
pressure drop. Figure 5 ma y therefore be used to evaluate 
the compressibility effects for other radiator pressure- 
drop definitions as well as for h 2 - P r > 7 



Once the ratio Ap/ipi is obtained from figure 5 

H 

for the design values of CDf. > Mr?> and n ^7n — * 

i c p6 rni r2 

the pressure drop for compressible flow may be obtained 
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from an evaluation of the pressure drop .for incompressible 
flow. For this purocse, use of tost data for the unheated 
radiator obtained with moderate rates of air flow (low 
Wach numbers) is generally desirable. This value of the 
pressure drop is then corrected to the design mass-flow 
and Reynolds number conditions in the usual manner without 
allowance for compressibility effects. With Api so 
determined, the correct pressure drop is computed from 

4 o = (^) C?: (7) 



COMPARISON 0 F THEORETICAL RESULTS WITH EXPERIMENTAL DATA 



A United amount of test data on the pressure drop 
In unheated tubes at Mach numbers up to the critical 
is available from references 3 and L. Similar data for 
heated tubes are not available. The available data will 
be analyzed and put into such form that direct comparison 
can be made with the theoretical results shown by the 
curves of figure 5( a )« 

Experimental data from reference 5 . - Reference J 
presents pressure"-^3rop data" for cylindrical tubes of 
circular section with rounded entries. nir was admitted 
from a large reservoir at atmospheric pressure. The 
pressure at the exit was progressively reduced until the 
limiting flow condition was attained. The measured mass- 
flow rate was expressed nondimensionally as m/m cr . This 
term is the ratio of the actual mass flow through the tube 
to the mass flow that would e;-:ist if sonic velocity were 
attained at the tube entrance. This parameter may be 
expressed in terms of the entrance Mach number as 



m 



1.7*28 M r2 



*cr 



Pressure data from reference 3 for 20-millimeter-diameter 
tubes with L/d = 10 and L/d = 60 are given in table I fe 
Values of M r2 corresponding to equally spaced values, of ^— , 

for which the pressure data were taken from reference 3, 
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are listed as column 2 of table I. The maximum values of 
m/m cr and M rp shown in table I were the highest that 

could be attained in the tests regardless of the outlet 
pressure and represent the condition at which sonic 
velocity was reached at the tube exits. 

In computing the values of Ap/&?. from the experi- 
mental data, equation 2(a) was used; that is 



and, since p 2 was a constant in the tests of reference 3, 



The experimental values of columns (1) and (J), were used 
to determine the value in column (k) of table I. The value 
of K was then chosen to make the experimental value 
of Ap/Ap^ at the lowest test speed (¥i r? = O.2.I4.Q) agree 



with the theoretical value obtained from figure 5(a) for 
this low speed. The resulting experimental values 
of Ap/Ap^ are shown in column 5. The theoretical values 

of Ap/Ap- shown in column 7 for comparison with these 
experimental results were taken from figure 5(a) for the 
values of M r ^ of column 2 and for values of C£ f . 

computed from figure 2 for 20-millimeter-diameter tubes 
with L/d = 10 and L/d = 60. Standard atmospheric 
entrance conditions and mass-flow values corresponding 
to m/m cr are assumed. The results shown in columns 5 
and 7 are plotted against M r2 in figure 6. 

The agreement between the pressure-droo ratios for 
the theory presented and for experiment will be seen to 
be unusually close. The limiting flow condition, in 
which sonic velocity was attained at the tube exits, was 
also accurately predicted by the theory. 
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Experimental data from reference L .- The data given 
in reference i| were obtained with a cylindrical pipe 
0.375 inches in diameter and ten feet long. Sxtremeiy 
high pressure drops were required to induce high-speed 
flow through the tube because of its unusually high L/d 
ratio of 320. Inlet pressures of several atmospheres 
were used. The tests consisted of measurements of 
pressure distribution along the pine for four combinations 
of inlet pressure and outlet pressure. The results of 
the tests were analyzed In reference I4. to determine values 
of the effective friction coefficient cf . The value 
of Cf was found to be Independent of the Mach number 
and to depend on only the Reynolds number, as assumed in 
the present report; furthermore, the values of Reynolds 
number and cf were virtually the same at all stations 
along the tube, as was assumed in the present analysis. 
Values of CBf^ from the test data agreed with 

theoretical values from figure Z to within I4. percent. 

In evaluating experimental values of the over-all 
presnure-droo ratio Ap/Ap^, Ap^ was calculated from 

equation (2a) as 




The basic data required for evaluation of Atu from 

this form of equation (2) are listed in reference 3. 
The over-all pressure drop Ap and the Mach number M r ^ 

were obtained directly from tabulations of the pressures 
and the Mach number in reference 3* ^ comparison of the 
test values of Ao/Apj_ with theoretical values taken 
from figure 5(a) is shown in the following table: 



i ■ ' ■" 


Experimental 






; Theoretical 


p 2 

(lb/sq ft) 


i 

j »r a 

1 
1 


i 


i 

(lb/sq ft) j 


Ar> 


I 
i 

i Ap 
-i 


7,1.22 

16,179 
17,607 


! 0.270 
1 .308 

! .32? 

i .297 


i s.8k 
1 k.9k 
\ 4-.13 
1 1+.12 


2,073 \ 
5 Ml 

11,891 : 
7,555 ! 


1.1+7 
s 2.02 

a 2.l8 

1.51 


| 1.51 
! a 2.09 

i a 2.11 

! i.l? 



a Sonic velocity attained at tube exit. 
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Good agreement between theory and experiment is shown. 
The accuracy of the theoretical curves of figure 5 has 
been verified for the conditions of extremely high inlet 
pressures, as well as for the more usual conditions 
(atmospheric inlet pressure) existing for the tests of 
reference 3, which are shown in table I and figure 6. 

A3 a matter of further interest, a check on the 
validity of the present theory may be made by computing 
the pressure distribution along the tube and comparing 
the results with the experimental data. These calculations 
were made by use of the present theoretical method, in 
which the pressure drops were computed for tubes of the 
same diameter (3/3 inch) as the test model but of lengths 
varying from 1 to 10 feet. The computed pressure drop 
for a tube 3 feet long, for example, was compared with the 
measured pressure drop from the tube entrance to the 
3-foot station. The theoretical pressures obtained from 
these computations are compared with experimental values 
in figure 7- Reasonably good agreement will be noted for 
all four test conditions. 

A principal assumption of the present theory was 
that the velocity profile across the tube section was 
uniform at all stations along the tube. The results of 
the comparisons in table I and in figures 6 and 7 indicate 
this assumption to be fully justified. In the analysis 
of reference J, in which account was taken of the 
velocity-profile shape, a similar conclusion was obtained 
theoretically. 



EXAMPLE OP USE 0? CURVES 



The use of the curves will be illustrated by 
calculating the pressure drop for an oil radiator for 
typical operating conditions. The dimensions of the 
radiator are as follows: 

Diameter, inches 

Duct diameter, inches 

Tube length, inches 

Tube diameter (inside), inches . . . 
Tube diameter (outside), inches . . 

Number of tubes 

Free-flow area, A r2 , square feet . 

Duct area, A2 and A*, square feet 

^CONFIDENTIAL 
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12 
12 
O.I96 
0.210 
1900 
0.397 

0.785 
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Basi c radiator performance dat a.- Performance data 
for the radiator are^ shown in figure 8 in the form that 
is most generally used. The air pressure-drop data given 
in curves of this type are frequently of limited value 
because the pressure drop is not adequately defined and 
because the test conditions under which it was measured 
are not stated. as previously discussed, the pressure 
drop may be measured by any one of at least four methods. 
It is obviously necessary to state which method was used. 
Pressure-drop data for the unheated radiator are valuable 
because they can be conveniently corrected to actual 
operating conditions by analytical means. Data for a 
specific heated condition are inconvenient to use because 
they must first be corrected to the unheated condition. 
The very wide variation of the heat-input parameter 
encountered in actual Installations makes it unlikely 
that the test values of this parameter will coincide with 
the required design value. 

Ass umed operating conditions .- The following assumed 
operaTTng conditions a op \y to a typical oil- cooler 
installation in an airplane flying in standard sea-level 
air at a speed of the order of 1+00 miles per hour: 

Air temperature at station 2, °P 88 

Air pressure at station Z, p£, 

pounds per square foot , absolute 2-55k 

Air density at station 2, pp, 

slugs per cubic foot 0.0C272 

Average oil temperature, °P <,.... 200 

Oil flow, pounds per minute 120 

Heat rejection, H, 3tu per second 75 

The heat rejected per minute per hundred degrees 
temperature difference between the oil and the inlet-air 
stagnation temperature is 

60 * 75 * 100 = ^ 02Q Btu per mlnute oer 100 o p 

200 — co 

With this value, the weight flow of air required is found 
in figure 8 to be ij.62 pounds oer minute. The mass flow 
is 

162 

m = 



60 x $2.2 



= O.239 slugs oer second 
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The required mass flow will deviate slightly from this 
value if the tube Reynolds number varies from the test 
conditions of figure 8. On the basis of the pro- 
portionality relationship between the heat-transfer and 
skin-friction coefficients, the change in skin-friction 
coefficient with Reynolds number (fig. 2) may be used to 
estimate the corresponding change in heat-transfer 
coefficient and hence to determine tho change in mass 
flow required. This effect is usually slight and, since 
only the mass-flow determination is involved, is considered 
to be outside the scope of the present report and will 
be neglected. 

Calculation of friction drag coefficient, C D .- 

_£i 

The tube Reynolds number is 
r = -JffiL- 



0.259 x 



0.19^ 
JL2_ 



0.397 * 3890 x 10- 10 



= 25,300 



From figure 2, for this value of Reynolds number and for 
L/d = 61.2, 

c f = O.OO57 

and 

C n = k x 61.2 x 0.0057 
f i 

1 = 1.1+0 

Calculation of Mach number .- The use of the curves 
of figure 5 requires the Much number at the tube 

entrance. Figure 3 gives K directly in terms of the 
known flow conditions ahead of the radiator; thus, 
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m?2 



0.5 



A r2 p 2 



0.239 x (jj.60 + 89) 
0.397 * 2^.0 ~ 

= 0.0554 



0.5 



With this value, the Mach number is found from figure 3 
to be 



M r2 = O.198 



Calcul at ion of heat-input paramet er.- The rate of 
heat input is specified on the curves oT figure 5 hy the 
nondimen clonal parameter H/c p gmT r 2 • This quantity is 

evaluated from equation (6) as 



H 



75 



c D 3mT r2 - °- 12 9? ^0.259 x 5^8 



0 



= O.O7I4. 



1 + 0.2(0.198) 



Compressibility effect .- The pressure-drop ratio 

Ao/Ap^ is taken from figure 5 for CDf. = 1 .i+0 , 

F 1 

Mrn = 0.1Q8, and ^ — = 0.07lj. # Interpolation between 

r 2 y c p gmT r2 

the curves of figure 5(b) and figure 5(c) is necessary 
because the heat-input parameter is intermediate between 
the values for these two parts of figure 5* The pressure- 
ratio value obtained is 



Ap 
APi 



- 1.17 



For this case, the effect of compressibility i s to 
increase the pressure drop 17 percent above the value 
for incompressible flow. In order to compute the actual 
pressure drop Ap, Ap^ must first be obtained. 



Calculation of AP4 



The basic incompressible-flow 



value of the pressure drop can be obtained from low-speed 
test data or from calculations based on the radiator- tube 
friction coefficients shown in figure 2. The method 
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involving the theoretical calculations requires the use 
of equation (2d) because the test ores sure -dron data for 
the present example is in terms of p^ - Pj • Prom 

Arp , 

equation (2d), therefore, since A 2 = A* and — ± = O.5O6, 

= 1 + 1.I4.0 - (0.506)* - 2^0.506 - (0,506) 

= 1.614- 

This commuted value of £0. /q.^ will now be compared with 
the value obtained from the test data of figure 8. 

For a pressure drop of 5 inches of water, the flow 
may safely be considered incompressible. At this test 
condition the weight flow is 197 pounds oer minute. 
Therefore, 

197 

m = 



32.2 x oO 
- 0.102 slugs per second 
and, for incompressible flow, 



1_ / 0.102 N 

2(0.002.576) \p.397 

13.9 



whence, for m = 0.102, 



'AoN _ 3 x 3-21 
~ 13.9 



= 1.87 
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This value is greater than the previously calculated 
value for m = 0.239 because of the lower tube Reynolds 
number at Ap = 5 inches of water (m = 0.102). The 
Reynolds number for this low mass flow is 



rad 

R = 7 



i r2 H 

0.196 
0.102 x — — 

0.397 * 3950 x 10" 10 
= 10,600 

The value of the tube friction-drag coefficient for this 
mass flow is, from figure 2. 



_ || x 12 x 0.0068 




= 1.66 

The value of CDf i for the design condition has 

been previously commuted as l.hQ. The correction to be 
aoolied to the low- speed test value of the pressure drop 
(see equation (2)) is q r2 (1.66 - 1.1+0); therefore, 




= 1.6l 

This result, based on the test data, compares very 
satisfactorily with the value 1.61+ previously computed 
from equation (2b) and equation (7)« 

The value of for the design condition is 

& Pi = 1.61+ (^r 2 ). 



CONFIDENTIAL 



MAC A ACR No. L^I23 



CONFIDENTIAL 



23 



\2 



Ap 



= 1.6k x 



V 0 _*Q7 / 



2 x 0,00272 \ 0.397/ 

= 110 Dounds per square foot 

Calculation of Ap .- The value of the pressure drop 
corrected for compressibility is, finally, 



Jp 
APi 



x 



= 1.17 x HO 

= 129 pounds per square foot 

Altitude effect. - The calculations rcade in the 
illustrative exa:.;nle for standard sea-level conditions 
have been roneattc. for standard conditions at an altitude 
of 30,000 feet. The san.e heat rejection, rate of oil 
flow, oil temperature, and airplane speed were assumed 
for both cases. The results obtained in the two cases 
are compared in the following table: 



^uanti ty 


Sea 


30,000 




level 


feet 


Air temperature, t2, °F 


88 


-19 


Air pressure, p£, lb/sq ft 


25^0 


795 


Air density, p£, slugs/cu ft 


0.00270 


0.00105 


Heat rejection, 3tu/min/l00°p 
Mass flow, slugs/sec 


1;020 


2055 


0.259 


0.070 


Reyno Id s number 


25,300 
1 .1+0 


8720 


Friction coefficient, CDf^ 


1.72 


Mach number, K r p 


O.I98 


O.165 


Heat-input factor, H/cpgmTr2 


0 . 07U 


0.314 


Pressure drop with compressibility 






neglected, lb/sq ft 


110 


29 


Pressure drop corrected for com- 






pressibility, Ap, lb/sq ft 


129 


ill 


Comoro 3 si bility correction, 






percent Apj_ 


17 
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The compressibility effect is shown to increase 
from 17 percent at sea level to i^l percent at 30,000 feet 
fnr'the typical oil-radiator installation assumed in the 
example. It is of interest to note that the oil-cooler 
tube" Mach number is somewhat less at altitude than at 
sea level. A very large increase in the heat-input^ 
parameter occurs, however, because of the decrease in 
mass flow required and the heating effect therefore 
becomes very large. 

Other tyces of cooler installation frequently require 
an increasing value of M r<2 with altitude. In the case 
of ethylene-glycol radiators, for example, the Mach number 
effect generally becomes very large at altitude and, in 
many cases, makes it impossible to obtain adequate 
cooling-air flow. 

Summary of method.- The steps required in the use 
of t he curves of figu re 5 to obtain the corrected value 
of the pressure drop are summarized as follows: 

(1) Determine the pressure, temperature, and density 
of the air ahead of the radiator and the heat to be 
dissipated. 

(3) Obtain the required mass flow of cooling air 
from the usual heat-dissipation data for the radiator. 

(3) Compute the Reynolds number of the tube flow 
and obtain the friction drag factor CD fi from figure 2. 

(il) Obtain the Mach number of the tube flow M r2 
from figure 5 . 

(p) Evaluate the heat-input factor H/cpgmTr 2 *>7 
equation ( 6 ) . 

(6) From figure 5 with these values of CDf ± , Mr 2 * 
and H/cpgmT r 2, obtain the pressure-drop ratio Ap/APi- 

(7) Compute APi from equations (2) or obtain it 
from low-speed test data corrected to the design 
Reyno 1 d s nuin be r . 

(8) Evaluate Ap by multiplying Api by the 
pressure-drop ratio Ap/Apj_. 
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CONCLUSIONS 



Simple curves are presented by which the basic 
nressure-drop characteristics of unheated tubular 
radiators can be corrected to operating conditions in 
which the radiator is heated and in which the Mach number 
of the tube flow is of appreciable magnitude. The 
following conclusions are indicated: 

1. The accuracy of the simple method presented for 
evaluating the compressibility effects is verified for 
unheated tubes for the complete range of attainable flow 
(tube exit Mach numbers from 0 to 1.0) by experimental 
data from two sources. 

2. The effects of compressibility (density change) 
on the pre 3 sure drop required to provide a given mass flow 
of cooling air through a tubular radiator are shown to 

be of appreciable magnitude under present-day operating 
conditions. Calculations for a typica 1 oil-cooler 
installation indicated that the compressibility effect 
increased the pressure drop by 17 percent at; sea level 
and by Ll percent at an altitude of 30,000 feet. 

Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va. 
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TABLE I 

ANALYSIS OF DATA FOR 20 -MILLIMETER -DIMETER 
TUBES FROM REFERENCE 3 



I Experimental 


l 

1 Theoretical 

* 


m 




Ap 


AE 
p 2 


An 


c Df . 


Ap 




°2 


(m/m cr ) 2 






Ap i 


L/d ~ 10 


0.1,0 

• *-h ) 

• 50 

.60 

.65 
.70 

• 75 
.80 

.85 
a -93 


0.2l|.0 
.273 

.307 
.3L2 
.;'7'J 

.la 8 

.1+58 
.502 

• 555 
.610 
.678 
.723 


0.050 
.065 

.o3o 
.100 
.120 
.145 
.175 
.210 

.255 

.315 
.400 

• 550 


0.312 
.320 
.320 
•33,1 

■M 

♦ 357 

• 373 

•m 
.L30 

.1x93 
.64° 


1.05 

1.03 
1.08 
1.11 
1.12 
1.16 
1.20 
1.26 

1.3I4. 
1-^7 
1.66 

a 2.l6 


0.171 
.166 
.16.1 
.153 

.154 
.151 
.148 
.11x6 

.11x3 
. 1I1 1 

.159 
.133 


1.05 
1.07 
1.09 
1.11 
1. 14 
1.17 
1.22 
1.27 
1.34 
1-45 

1.65 
2.01 


L/d = 60 


0 .1+0 

•45 
.50 

.60 

.65 

.70 

.75 
a .77 


O.240 

.273 
.507 

• 3k2 

.578 
.£18 

.458 
.502 
.522 


0.085 
.110 

.li|0 

.175 

.220 
.275 

.L60 
.630 


0.530 

• 543 
.560 

• 578 
.612 
.650 

!8l^ 
1.062 

« 


1.10 
1.13 

1.17 ! 

1.20 

1.27 

1.35 

l.fo 
1.70 
a 2.22 


1.026 
.996 
.966 
.948 
.92ix 
.906 
.888 
.876 
.870 


.. 1.10 
1.15 
1.17 
1.22 
1.28 
1.36" 

1.5-7 
1-73 
a 2.20 



a Sonic velocity attained at tube exit. 
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Fig. 2 NACA ACR No. L4I23 




NACA ACR No. L4I23 



Fig. 
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Figure 6 . — Van at /on of static- pressure - drop ratio 
with Moch number. 
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Figure 5 . - Continued. 
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Figure 5 r Concluded. 



Fig. 6 
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Fiqure 6 . ~ Comparison of fheoref/cal variation of 
pressure -drop raf/o w/fh experimental data 
obtained from reference 3. 
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